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DEVELOPMENT OF A COMPUTER PROGRAM TO OBTAIN ORDINATES 
FOR NACA 4-DIGIT, 4-DlGIT MODIFIED, 5-DIGlT, 

AND 16-SERIES AIRFOILS 

Charles L. Ladson and Cuyler W. Brooks, Jr. 

Langley Research Center 

SUMMARY 

The analytical design equations for both symmetrical and cambered airfoils in the 
NACA 4-digit, 4-digit modified, 5-digit, and 16-series airfoil families have been reviewed. 
A computer program has been developed to calculate rapidly the ordinates and surface 
slope for these airfoils and the program is included as an appendix to this report. Pro- 
visions are made in the program to combine basic airfoil shapes and camber lines from 
different series so that nonstandard airfoils can also be generated. The program also 
produces plots of the nondimensional airfoil ordinates and a punch card output of the ordi- 
nates in the inp It format of a readily available program for determining the pressure dis- 
tributions of arbitrary airfoils in subsonic potential viscous flow. 

INTRODUCTION 

During the 1930’s several families of airfoils and camber lines, all of which have 
analytic expressions for the ordinates, were developed by the National Advisory Committee 
for Aeronautics (NACA). These include the NACA 4-digit airfoils (ref. 1), 4-digit modified 
airfoils (ref. 2), 5-diglt airfoils (ref. 3), and 16-series airfoils (refs. 4 and 5). Many of 
these airfoil shapes have been successfully used over the years as wing and tall sections 
for general aviation as well as military aircraft. Some have been and are still being used 
as sections for propellers and helicopter rotors. 

Numerous specific airfoils of these series have been computed and ordinates pub- 
lished over the year?:. However, when performing parametric studies on effects of such 
variables as thickness, location of maximum thickness, leading-edge radius, amount and 
location of maximum camber and others, it is not always easy to obtain the ordinates of 
the desired shapes rapidly. Because these airfoils all have analytic solutions for the 
ordinates, both with and without camber, a computer program can be written to provide 
the exact ordinates rapidly and at a low cost. An attempt to do this was made in refer- 
ence 6, but some limiting assumptions were made so that exact results are not provided 

for some airfoils. 



The purpose of this paper is to review the design parameters for all these airfoils 
and to describe a computer program which will generate exact ordinates for all airfoils 
of these series with an acceptable expenditure of computer time. The program will also 
allow combination of any airfoil and any camber line so that many nonstandard airfoils can 
be described. 


SYMBOLS 

When two symbols are given for a concept, the symbol in parenthesis is that used in 
the computer program and on computer- generated plots, 

A camber line designation, fraction of chord from leading edge over which 

design load is uniform 

aQ,aj,a2,a3,a^ constants in airfoil equation 

bQ,bj^,b2 constants in camber line equation 

c (C) airfoil chord 

(Cl) (CLI) design section lift coefficient 

V 'design 

dQ,dj,d2,d3 constants in airfoil equation 

I leading-edge index number 

constants 

m chordwise location for maximum ordinate of airfoil or camber line 

p maximum ordinate of 2-digit camber line 

R radius of curvature 

r chordwise location for zero value of second derivative of 3- digit camber 

line equation 

t thickness 


X (X) 
2 


distance along chord 


y (Y) airfoil ordinate normal to chord, positive above chord 

5 local inclination of camber line 

Subscripts: 

cam cambered 

I (L) lower surface 

te leading edge 

N forward portion of camber line 

T aft portion of camber line 

t thickness 

u (U) upper surface 


ANALYSIS 


The design equations for the analytic NACA airfoils and camber lines have been pre- 
sented in references 1 to 5. They are repeated herein to provide a better understanding 
of the computer program and indicate the use of different design variables. A summary 
of some of the design equations and ordinates for many airfoils from these families is 
also presented in references 7 to 9. 


Thickness Distribution Equations 

4-digit .- Ordinates for the NACA 4-digit airfoil family (ref. 1) are described by an 
equation of the form: 



The constants in the equation were determined from the following constraints: 
(1) Maximum ordinate; 


0.30 
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(2) Ordinate at trailing edge: 

1 = 1.0 I = 0.002 

(3) Trailing- edge angle; 

2 = 1.0 g = 0.234 

(4) Nose shape: 

1=0.1 1=0.078 

The coefficients listed below were determined to meet these constraints very closely; 
ag = 0.2969 

aj = -0.1260 

a2 = -0.3516 
a^ = 0.2843 
a4 = -0.1015 

To obtain ordinates for other thickness airfoils in the family, the ordinate .5 for the 
0.20-thickness-ratio model are mtiltiplied by the ratio (t/c)/0.20. The leading-edge 
radius of this family is defined as the radius of curvature of the basic equation evaluated 
at § ■ 0. Because of the term aQ^^x/c in the equation, the radius of curvature is finite 
at this point and can be shown to be ao^/2. Thus, the leading-edge radius varies as 
the square of the airfoil thickness-chord ratio because the thickness varies linearly with 
the a constants. To define an airfoil in this family, the only input necessary to the com- 
puter program is the desired thickness-chord ratio. Symmetric airfoils in this family are 
designated by a 4-digit number, that is, NACA 0012. The first two digits indicate a sym- 
metric airfoil and the second two, the thickness-chord ratio. 

4 -digit modified .- The design equation for the 4-digit airfoil family was modified 
(ref. 2) so that the same basic shape was retained but variations in leading-edge radius 
and chordwise location of maximum thickness could be made. Ordinates for these airfoils 
are determined from the following equations: 


4 


I 


1 


From leading edge to maximum thickness, 

4 =“o#+m(|)+» 2 (I) +“3(1) 

Prom maximum thickness to trailing edge, 

±1 = do + dj^l - + d2^1 - I) +^3^1“ 

The constants in these equations can be determined from the following constraints: 


(1) Maximum ordinate: 




1 - 0.1 


0 

dx 


(2) Leading- edge radius: 


0 


R = 


ao 


c - - 2 

(3) Radius of curvature at maximum thickness: 


X 


= m 


R = 


(1 - m)^ 


2di(l - m) - 0.588 


(4) Ordinate at trailing edge: 


1=1.0 


i = do = 0.002 


(5) Trailing- edge angle: 


1= l.C 


^ = di = f(m) 
dx ^ 


Thus, the maximum ordinate, slope, and radius of curvature of the two portions of the 
surface match at f “ m. The values of d^ were chosen, as stated in reference 2, to 
avoid reversals of curvature and are given in the following table; 


m 

dl 

0.2 

0.200 

.3 

.234 

.4 

.315 

.5 

.465 

.6 

.700 
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By use of these constraints, equations were written for each of the constants ^except 
dQ and d^) In the equation for the airfoil family and arc included in the computer pro- 
gram. As in the 4-diglt airfoil family, ordinates vary linearly with variations in thickness- 
chord ratio and any desired thickness shape can be obtained by scaling the ordinates by the 
ratio of the desired thickness ratio to the design thickness ratio. 

These .airfoils are designated by a 4-digit number followed by a dash and a 2-digit 
number (that is, NACA 0012-63). The first two digits are zero for a symmetrical airfoil 
and the second two digits indicate the thickness -chord ratio. The first digit after the dash 
is a ’eading-edge radius index number, and the second is the location of maximum thick- 
ness in tenths of chord aft of the leading edge. The leading-edge index is an arbitrary 
number assigned to the leading- edge radius in reference 2 and is proportional to Uq. An 
index of 0 indicates a sharp leading edge (radius of zero) and an index of 6 corresponds to 
aQ = 0.2969, the normal design value for the 4-digit airfoil. A value of leading-edge index 
of 9 for a three times normal leading-edge radius was arbitrarily assigned in reference 2. 
Values of leading-edge radius for various values of the index number and thickness-chord 
ratio are listed in table I and plotted in figure 1. The computer program is written so that 
the desired value of leading-edge radius is the input parameter. The value of a^ is then 
computed in the program. The index number is only used in the airfoil designation. 

16-series .- The NACA 16-series airfoil family is described in references 4 and 5. 
Although not directly stated in the references, it will be noted from the equation for the 
ordinates in reference 5 that this series is a special case of the 4-digit modified family. 
The 16-series are thus defined as having a leading-edge index of 4 and a location of max- 
imum thickness at 0.50 chord. The designation NACA 16-012 airfoil is equivalent to an 
NACA 0012-45. The computer program does not have separate inputs for the 16-series so 
that the 4-digit modified series must be used to obtain ordinates for these airfoils. 

Camber- Line Equations 

2-digit .- The NACA 2-digit camber line is described in reference 1. This camber 
line is formed by two parabolic segments which have a general equation of the form 

^ = bo + The constants for the two equations are determined from the fol- 

lowing boundary equation -j: 

(1) Camber-line extremities: 

c c 

1=1.0 |.o 
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(2) Maximum ordinate; 



X 

c 


m 


1 = 
c 


P 


dx 


= 0 


From these conditions, the camber-line equations then become 



forward of maximum ordinate and 


1 = 
c 


(1 - m)^ 


(1 - 2m) + 2m(|) 


aft of the maximum ordinate. Both the ordinate and slope of the two parabolic segments 
match at ^ = m. This camber line is designated by a two-digit number and, when used 
with a 4-digit airfoil, would have the form NACA pmXX where p is the maximum camber 
in percent chord; m is the chordwlse location of maximum camber; and XX is the 
airfoil thickness in percent chord. Tables of ordinates for some of these camber lines 
are tabulated in references 8 and 9. The ordinates are linear with amount of camber and 
these can be scaled up or down as desired. 


3-digit.- To provide a camber line with a very far forward location of the maximum 
camber, the 3-diglt camber line was developed and presented in reference 3. This camber 
line is also made up of two equations so that the second derivative decreases to zero at a 
point r aft of the maximum ordinate and remains zero from this point to the trailing 
edge. The equations for these conditions are as follows: 

From 1=0 to | = r, 



From ^ = r to ^ = 1.0, 
c ^ 



0 
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The design criteria are as follows: 
(1) Camber-line extremities: 




s-l.o 

(2) At junction point: 



The equation for the camber line then becomes 



from 1=0 to 1= r and 



from 3 * r to ^ = 1.0. These equations were then solved for values of r which would 
give longitudinal locations of the maximum ordinate of 5, 10, 15, 20, and 25 percent chord. 
The value of was adjusted so that a theoretical design lift coefficient of 0.3 was 
obtained at the ideal angle of attack. The value of kj^ can be linearly scaled to give any 
desired design lift coefficient. Values of k^ and r and the camber-line designation 
were taken from reference 3 and are presented in the following table: 


Camber-line 

designation 

x/c for maximum 
camber, m 

r 

^1 

210 

0.05 

0.0580 

361.400 

220 

.10 

.1260 

51.640 

230 

.15 

.2025 

15.957 

240 

.20 

.2900 

6.643 

250 

.25 

.3910 

3.230 


The first digit of the 3-digit camber-line designation is defined as two-thirds of the design 
lift coefficient, the second digit as twice the longitudinal location of maximum thickness 
in tenths of chord, and the third digit of zero indicates a nonreflexed trailing edge. 
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3-digit reflex .- For some applications, for example, rotorcraft main rotors, it may 
be desirable to produce an airfoil with a quarter-chord pitching -moment coefficient of zero. 
The three-digit reflexed camber line was thus designed to have a theoretical zero pitching 
moment as described in reference 3. The forward part of the camber line is identical to 
the 3-digit camber line but the aft portion was changed from a zero curvature segment to 
a segment with curvature. The equation for the aft portion of the camber line is expressed 

by s k 2 (“ - r^ By using the same boundary conditions as were used for the 3-digit 

earner line, the equations for the ordinates are 


y . 1 

■ 6 






from 1=0 to I = r and 


c 6 1 




for - = r to - = 1.0. The ratio kg/kj is expressed as 


^ 3(r - m)^ - r^ 

kj “ 1 - r 


Values of kj, ^2/^V ^ several camber-line designations from reference 2 are 

presented in the following table: 


Camber-line 

designation 

x/c for maximum 
camber, m 

r 

^1 

HlH 

221 

0.10 

0.1300 

51.99 

0.000764 

231 

.15 

.2170 

15.793 

.00677 

241 

.20 

.3180 

6.520 

.0303 

251 

.25 

.4410 

3.191 

.1355 


The camber-line designation for this camber line is identical to that for the 3-digi.t camber 
line except that the last digit is changed from 0 to 1 to indicate the reflex characteristic. 

6- and 6A-series .- The equations for the 6-series camber lines are presented in 
reference 8. These camber lines are a function of the design lift coefficient 
and the chordwise extent of uniform loading A. These 16-series cambered airfoils (ref. 4) 
are derived by using the A = 1.0 camber line of the series. These equations have been 
programed for use with 6-series airfoils in reference 10 and that part of the program has 
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been Incorporated Into the present program. As was the case in reference 10, the pro- 
gram Is capable of combining up to 10 camber lines of this series to provide many types 
of loading. 


Calculation of Cambered Airfoils 

To calculate ordinates for a cambered airfoil, the desired mean line is first com- 
puted and then the ordinates of the symmetrical airfoil are measured normal to the mean 
line at tlie same chord station. This procedure leads to a set of parametric equations 
where (y/c)^, (y/c)(,am> ^ are all functions of the original independent variable 

x/c. The ordinates on the cambered airfoil (x/c)^ and (y/c)^ are given by 



where 6 is the local inclination of the camber line and (y/c)^ is assumed to be nega- 
tive to obtain the lower surface ordinates (x/ c)^ and (y/ c)^. This procedure is also 
described in reference 1. The local slopes of the cambered airfoil can be shown to be 



and 



by parametric differentiation of (x/c)^^^ and (y/c)^^^ with recpect to the original x/c 
and use of the relationship 

/dy\ d(y/c)u/d(x/c)u 
\^u ' d(x/c)u/d(x/c) 

Although specific camber lines are generally used with specific thickness distri- 
butions, this program has been written in a general format. As a result, any camber line 
can be used wdth either type thickness distribution so that any shape desired can be 
generated. 
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RESULTS AND DISCUSSION 


I 

1 


Program Capabilities 

The computer program which was developed to provide the airfoil shapes described 
by the equations in the analysis section is listed in the appendix. The output of the pro- 
gram consists of tabulated ordinates, computer-generated plots of the non dimensional 
ordinates, and punched card listings of the ordinates. The punched cards are in the format 
of the input of the program described in reference 11 so that pressure distributions over 
the generated shape may be readily obtained. To show graphically the capabilities of the 
program, sample computer plots of several airfoil shapes are presented in figures 2 to 10. 

Figures 2 and 3 illustrate possible variations in the 4-digit airfoil family, figure 2 
showing variations in thickness -chord ratio for symmetrical airfoils and figure 3 showing 
variations in the amount of camber for a fixed thickness-chord ratio and location of maxi- 
mum camber. Figures 4 and 5 illustrate possible variations in the 5-digit airfoil family. 
Variations in the longitudinal location of maximum camber are shown in figure 4 and a 
comparison of the same airfoil with nonreflex and refle>: camber lines is shown in figure 5. 
Examples of the 4-digit modified-series are shown in figure 6 for symmetrical airfoils and 
in figure 7 for cambered airfoils. The symmetrical airfoils have variations in the longi- 
tudinal position of maximum thickness whereas the cambered airfoils show variations in 
the longitudinal position of maximum camber. 

Examples of 16-series airfoils (which, as previously noted, are special cases of 
4-digit modified airfoils) are shown for symmetrical and cambered sections in figures 8 
and 9, respectively. Figure 10 presents an example of a combination of a 4-digit modified 
airfoil with a combination of two 6-series camber lines to give an aft-loaded section. This 
is shown to give an indication of the types of sections which may be generated by combi- 
nations of various thickness distributions and types of camber lines. If a thickness -chord 
ratio of 0.0 is specified in the input to the program, the shape of just the camber line or 
combination of camber lines is computed. The results of this procedure are shown in 
figures 11 and 12. 


Sample Output Tabulations 

Sample computed ordinates for both a symmetric and a cambered airfoil are pre- 
sented in tables II and III, respectively. Printed at the top of the first page for each table 
is the airfoil and camber-line family selected, the airfoil designation, and a list of the 
input parameters for both airfoil shape and camber line. lor the 4-digit modified airfoil 
family, the coefficients of the airfoil equation are also listed for a shape with a thickness- 
chord ratio of 0.20. Both nondimensional and dimensional ordinates are listed. The 
dimensional quantities have the same units as the input value of the chord, which is also 
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listed at the top of the page. First and second derivatives of the surface ordinates are 
also presented for symmetrical airfoils, but only first derivatives are tabulated for tlie 
cambered airfoils. 


Accuracy of Results 

All the airfoils and camber lines generateo by this program are defined by closed 
analytical expressions and no approximations have been made in the program. Thus, all 
results are exact. Many cases have been run and compared with previously published 
results to check the procedure and in all cases the comparisons were exact except for 
occasional differences in the last digit due to rounding differences. 

CONCLUDING REMARKS 

The analytic design equations for both symmetrical and cambered airfoils in the 
NACA 4-digit, 4-digit modified, 5-digit, and 16-series airfoil families have been reviewed. 
A computer program has been developed to calculate rapidly the ordinates and surface 
slope for these airfoils and the program is included as an appendix to this report. Pro- 
visions are made in the program to combine basic airfoil shapes and camber lines from 
different series so that nonstandard airfoils can also be generated. The program will also 
produce plots of the nondimensional airfoil ordinates and a punch card output of the ordi- 
nates in the input format of a readily available program for determining the pressure dis- 
tributions of arbitrary airfoils in subsonic potential viscous flow. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
August 29, 1975 
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APPENDIX 


COMPUTER PROGRAM FOR ORDINATES OF ANALYTICAL NACA AIRFOILS 

The program presented herein is written in the Langley Research Center version of 
FORTRAN IV and has been used on the Control Data series 6000 computer systems. Both 
the computational program and a plotting program are presented, although the plotting rou 
tine is included as a guide for users only. Several unlisted subroutines are used in the 
plotting program. The computational program requires about 46000g storage locations, 
and requires about 8 seconds to compile and about 1.5 seconds to execute each case on the 
Control Data 6600 computer system. 

Card Input Format 

The input to the program is in a card format as follows: 

CARD 1 - Number of ordinates to be output on punched cards: (Maximum of 32) (right 
justified in columns 1 to 3). 

CARDS 2,3,4, and 5 - Chordwise location of ordinates to be output on punched cards. 

(Columns 1 to 10, 11 to 20, etc., with decimai point.) 

CARD 6 - Tabulated data printout airfoil title card. Any designation may be used in col- 
umns 1 to 80. 

CARD 7 - Airfoil thickness series and camber-line series designations are as follows: 


NACA airfoil family 

Card designation* 

Columns 

4 -digit 

4 -digit modified 

4 -DIGIT 
4-DlGITMOD 

1 to 7 
1 to 10 


Camber line 

Card designation* 

Columns 

NACA 2-digit 

2 -DIGIT 

11 to 17 

NACA 3-digit 

3 -DIGIT 

11 to 17 

NACA 3-digit reflex 

3-DIGITREF 

11 to 20 

NACA 6-series 

6 -SERIES 

11 to 18 

NACA 6A-series 

6A-SERIES 

11 to 19 


*These are hollerith cards; designations must be in 


exact columns. 
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CARD 8 - Airfoil thickness distribution parameter card. (Note that cards 3 to 7 are in 
floating-point mode. Numbers are entered with a decimal point.) 


Description 

Variable 

Columns 

Thickness -chord ratio of airfoil 
(i.e., 0.120) 

TOC 

1 to 10 

Leading-edge radius to chord ratio. 
Not used with 4-digit but must be 
used with 4 -digit modified 

LER 

11 to 20 

Basic chordwise increment in x/ c 
for computing ordinates. Usually 
set to 0.01 

DX 

21 to 30 

Model chord used for listii^ ordinates 
in dimensional units 

CHD 

31 to 40 

Nondimensional chordwise location 
of maximum thickness. Used for 
4 -digit modified airfoils only 

XM 

41 to 50 

Trailing-edge slope of 4-digit modified 
airfoils. Take values from text or 
reference 2 or input 0.0 and approxi- 
mate value from equation in reference 7 
will be used 

D1 

51 to 60 


CARD 9 - Airfoil camber-line parameter card. Set all values equal to 0.0 for a symmet 
rical airfoil. 


Camber line 

Description 

Variable 

Columns 

2 -digit 

Maximum camber ordinate to chord 
ratio (i.e.. 0.04), p 

CMB 

1 to 10 


Longitudinal location of maximum 
camber position (i.e., 0.40), m 

CM 

11 to 20 

— ——I 

3 -digit 

Value of k^ from text or reference 3 
which varies linearly with design 
lift coefficient 

CMB 

1 to 10 


Value of r from text or reference 3 
which is a function of the longitudinal 
location of maximum camber 

CM 

11 to 20 
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Camber line 
3 -digit reflex 


6 series and 
6A-series 


APPENDIX 

Description 

Value of from text or reference 3 
for reflex airfoils which varies 
linearly with design lift coefficient 
Value of r from reference 3 for 
reflex airfoils which is a function 
of the longitudinal location of 
maximum camber 

Value of kg/kj from reference 3 for 
reflex airfoils which is a fvmction of 
loi^itudinal location of maximum 
camber 

Design lift coefficient (i.e., 0.20) 
Camber line chordwise loading 
(use 0.80 for 6A- series) 

Number of camber line to be 
summed (if only one, leave 
blank or insert 1.0) 


Variable Columns 


CMB 


K20K1 


1 to 10 


11 to 20 


21 to 30 


CLl 1 to 10 

A 11 to 20 

CMBNMR 21 to 30 


CARDS 10, 11, and 12 — Up to nine additional camber lines may be summed on these cards 

for the 6-series camber line. These cards are not necessary for 
only one camber line. 


Camber line 

Description 

Variable 

Colmnns 

6-series 

Design lift for second camber line 

CLI 

1 to 10 


Loadii^ for second camber line 

A 

11 to 20 


Design lift for third camber line 

CLI 

21 to 30 


Loading for third camber line 

A 

31 to 40 


Design lift for fourth camber line 

CLI 

41 to 50 


Loading for fourth camber line 

A 

51 to 60 


Design lift for fifth camber line 

CLI 

61 to 70 


Loading for fifth camber line 

A 

71 to 80 


CARD 13 - Title card for use in plot of airfoil ordinates. Any designation may be used In 
columns 1 to 80. 
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Program Listing 

A program listing follows: 




PROGRAM I N ( I NPut ^OUTPUT t TAPES* I NPUT • T APE6*0UTPUT tPUNCH ) 

DIMENSION XU(aOO), XL (200)* YU(200)* YL ( 200 > 

COMMON /MAIN/ YST ART ( 3 ) • CHD * KON 
DIMENSION C0EFFS(4) 

DIMENSION XA<32)* XAU(32)* YAU<32)* XAL ( 32 ) * YAL(32)* NAME<2) 
DIMENSION CLI(IO)* A(10)* TANTHO(IO)* YCMB(IO)* TANTHdO)* YCP2 < 1 0 
I ) * IF6XA( 10) 

INTEGERS TITLE (8) 
complex roots (3)* TEMP (8) 
equivalence (CLI ( I ) fCMB) 
integer SERIEC 
PEAL K20KI 

integer PROFILE* camber 

YSTART( 1 )=1 #0 

YSTART(2)*4#0 

YSTART(3)s7«0 

K4D=10h4-DIGIT 

K4DMO0* 10H4-DIGITMOD 

K2DalOM2-OIGIT 

K3DB10H3-DIGIT 

K3DREFS10M3-DIGITREF 

K6S»10H6-SERIES 

K6ASslOH6A-SERIES 

KONsO 

input parameters normalized by the chord (CHD) 

TOC - T/C* THICKNESS* RLE - LEADING EDGE RADIUS* XM - X(YMAx)/CH0R 

DX - INTERVAL/CHORD* CHD - CHORD IN DESIRED UNITS 

CMB - CAMBER CONSTANT Kl* CM - X(MAX CAMBER ) /CHORD 

CALL PSEUDO 

call LEROY 

read (5«S90) N* (XA ( I ) • I * 1 *N ) 

01 *0*0 

20 READ (5«600) ( T I TLE ( I) * I = 1 * 6 ) 

IF (ENPFILE 5) 30*40 
30 call CALPLT (0*0*999) 

STOP 

40 CONTINUE 

read (S*600) profile*camber 

KON*KON«M 
1 CKYa0*0 
FRAC=l .0 
Dial *0 

PRINT 62<?* PROFlLE*CAMBER 
PRINT 610* (TITLE ( I ) * I *1 *8) 


A 30 
A 40 
A 50 
A 60 
A 70 
A 80 
A 90 
A 100 
A 110 
A 120 
A 130 
A 140 
A ISO 
A 160 
A 170 
A 160 
A 190 
A 200 
A 210 
A 220 
A 230 
A 240 
A 250 
A 260 
A 270 
A 260 
A 290 
A 300 
A 310 
A 320 
A 330 
A 340 
A 350 
A 360 
A 370 
A 380 
A 390 
A 400 
A 410 
A 420 
A 430 
A 440 


IK (PWOKIl F, to. lnM4. DIGIT ) READ (b*S3(iJ TOC*RLF *DX*CwD»R*H*H*NA A 4 SO 

IMt(l) A 

XM = O.S ^ 

IK (PPOKII.F.E0.l(^H4*0lblTMOD) READ (b*^TU) ToC * RLF • DX * CHD • XM ♦ D l , H * A 480 

iNA'Vlt (1) ^ 

It- (CAMbE;^.EQ*iaH/-nlOlT .0^? •CAMrttp . Fo* 1 T ) rE AD A SQO 

10) CMR*CM*PB»HB*HB*RH*riH-,N6ML (/) A MU 

IK (Cr^riFP.EOdOHT-nlOlTRFF > (^>t630) CMB * CM * K 20K 1 * HB • RP ♦ 8^ • HR A B?U 

1 fNAMF (?) ^ 

U (CA'^HE^.EQ* 10hb-«;FP1FS .op . C AmHER , FO • 1 OHbA-btP IFS ) REAO (S.63 A R4U 

10) Cl I (1) *(i(l) •CMHNMR*bH,WB*dH*8H*NAMF(2) A RSU 

IF (CA‘1rtt*^.EQ*lOH^-SEHlFS .OP * C AHHtH . FO . 1 OHf>A-St R T KS ) 1 Cky*CMPNm A S6U 

I A B 7 0 

IK (IC^'r.DT.l) (b*640) ( CL I ( I ) f A ( D * 1 =2 , I CK Y ) AS8U 

It- (CtMbt^-'.FOtlOHB-BEPlKS .OW.CAMHLK.FQ. InHbA-St PIKS ) C»"R = CLT ( I ) A ^90 

PKlMT H-iO. NlA^h A AOO 


I 


i 
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If- (ICkY.LE.D ICKY = 1 

t = , 

UO SO I =l t 10 

lf-Sl(A< I )=o 

SO CONTiNUt _ . . ^ 

II- (PPOf- tl.F .to. I 0H^*-.Dll»l T ) point t OCf OLE •OX tCMO 

If- (PPOF II.F.tO. IOH 4 -.OI 0 I TMOlJ) point b70» TOCfOl.Ef(>X •CMOtXM,Ol 
If- (CA'^ttEO.EO* 10H?-nlOlT ) POINT OHO, C^'iR.CM 

U (C^^^rtE^v.EO. lOHl-niGl T ) POINT O'VO# C^^'RiCM 

Ih (C^'^oEP.tO* lOH'T-niOlTOt.F ) PPINT ^00* CMR • CM • K jj-OK I 

If- <CAMrfEo.EO* lOH^-qEOltS .oO»CAMbLW.E0* lOHbA-StPlFS ) PWImT 7lQ 

IF (CA^^bEO.EO. 10 H(S-c;E01Es .U« •CAMBlP .FO . 1 0H6A-st ^T Es > P0I^»T 7?0^ 

1 (a.i(i)fMnti = i»icKY) 

IT (TOC.LT.F.)POOFILEsK4U 

It (PROFILE. to. lOHA-OlbiT ) ^0 TO 70 

computed constants 

AU=SOOT (2.0ttRLE) oO.?/TOC 
00 = 0 . 00 ^? 

It (01 .GT.0.0) on in bO 

Ol = , (<?. 24 - 5 . 4 . !?• 1 1 .-0,878^ XM) 

bO CONTINUE 

U3=(3.<>01-0.SH8/(l.-XM) )/(3.*(l.-XM)o*?) 

02=-l .S« n .-XM) «n3-,5«Ul/ ( 1 .-XM) . ^ , w. 

AJ = 0« 1 /XMo«3* (^.^^01^^ ( 1 .-XM) -O.'iHH) / (2.»XM<=- ( 1 . -XM ) «»? ) - 3 . * A 0 / ( 8 . » X ^ 

l<K»a.5) 

A«r=-0. l0/XMO«2».S<>An/XM<»*l .5-?.«XM«A3 

A1 =-,5»A0/XM#*#5-2**XM*A2-3«*XM**2#A3 

RC IS RADIUS OF CURVATURE AT X=*XM 

RC=( ( 1 •-XM)**2/ <2**01*< I •-XM)-0«588) )*«2/tOC 

PRINT 730* AO *A1 *A2*A3*00*D1 *02*03* RC 

PROFILE* X LE XM 

70 continue 

IF (ABS(CM3)#LE*0*l*»6) PRINT 740 

IF (AeS<CM0)*GT.O*l»*6) PRINT 750 

X = 0«0 

Y = 0#0 

XC*0#0 

YC=0*0 

XU< I )«0.0 

YU( 1 )tt0«0 

XL (1 ) * 0 • 0 

YL(1 )*0#0 

xuc*o.o 

YUC«0»0 

XLC*0«0 

YLC*0«0 

XAU( 1 )*0«0 

YAU( 1 )aO«0 

XAL( 1 )=0«0 

YAL< I )aO*0 

K = 2 


IF 

(CAMBet)*EQ* 10H2-DIGIT 

) 

GO 

TO 

8C 

IF 

«CAMBER.e:Q.l«H3-0ICIT 

) 

GO 

TO 

90 

IF 

<CAHB£R*E0* tUH3-0iGITREF) 

GO 

TO 

100 

IF 

(CAMBER. EQ» 10H6-SERIES 

) 

GO 

TO 

110 

IF 

(CAMBER.EQ. JOHSA-SERIES 

) 

GO 

TO 

1 10 


PRINT 760 
GO TO 190 


A SIP 
A SPU 
A S3-J 
A S40 
A SSP 
A SSO 
A S70 
A S80 
A S90 
A 70u 
A 710 
A 720 
A 730 
A 74u 
A 7Al 
A 7SO 
A 760 
A 770 
A 78U 
A 790 
A AOO 
A AlO 
A A20 
A R30 
A A4 0 
A A50 
A A60 
A A70 
A A80 
A 890 
A 900 
A 910 
A 920 
A 930 
A 940 
A 950 
A 960 
A 970 
A 960 
A 990 
A 1000 
AlOlO 
A1020 
A 1030 
A 1040 
AiO*SO 
A 1060 
A 1070 
AlOBO 
A 1090 
A1 100 
AlllO 
A1 120 
A1 130 
Al 140 
A1 150 
A1160 
Al 170 
Aiiao 
Al 390 
A1200 
A1210 
A 1 220 
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eO tANTHO«2**CM0/CM A1230 

IF (A0S(CMa>«LT.E) TANTHOaE A1240 

YP»lO«#*lO A1250 

YPP«lO.*#lO A1260 

YUPe-l /TANTHO A 1270 

YLPs- 1 /TANTHO A 1 200 

GO TO 190 A1290 

9C TANTH0sCMB*CM#»2<M3«0-CM)/6*0 AI300 

IF ( ABS(CMB)«LT*E) TANTHO=E Al3lO 

YPslO.<HHO A1320 

YPPslO»**lO A1330 

YUPa- I /TANTHO A I 340 

YLPa- I / T ANtHO A 1 350 

GO TO 190 A I 360 

100 TANtHOsCM8*(3«*CM**2-K20Kl* ( I ~CM ) 3 ) /6 A 1 370 

IF (ABS(CMa)*Lt*E) TANTHOsE A 1 380 

YP=lO«#»lO A1390 

YPPalO«**lO A1400 

YUP*“1/TANTH0 A141 0 

YLPs- 1 /TANTHO A 1420 

GO TO 190 A1430 

I 10 L*0 A1440 

CLISaCLIdl A 1450 

AS=A ( I ) A1460 

120 LsL+l A1470 

A ( I ) sA (L) A 1480 

CLl ( 1 )=CLI (L) A1490 

Ks2 A I 500 

U=0#005 A151 0 

V«- ( A-U )/ABS< A-U ) A1520 

OMXL*(l t-Ul^ALOGI 1 .-U) A1530 

AMXL*(A-U)*ALOG( ABS(A-U) ) A 154 0 

0MXL1*-AL0G( 1 •-U)-l • A 1550 

AMXL13-AL0G(ABS<A-U) ) + V a 1560 

0MXL2*l •/< 1 *-U) A 1570 

AMxL2s-V/ABS(A-U) A 1580 

IF (A.LT»E#OR«ABSn •-AI.LT^E) GO TO 130 A1590 

Gs-( A*»2* ( #5 <^AlOG (A )-0*25 ) + 0«25) / ( 1 •-A ) A I 600 

Qa I , 0 A1 61 0 

H« (0*5» ( 1 #-A )4*2*AUOG ( I »-A )-0»25# ( 1 •*-A )#*2 ) / ( I »-A )+G A I 620 

Z=*5» ( A-U )»AMXL-*5* ( 1 • -U ) ♦OMXL- • 25* < A-U ) **2+ • 25* ( I •“U)**2 A 1630 

21 a«5* ( (A-U)*AHXL1 -AMXL- ( 1 • -U ) *OMXL 1 +OMXL+ ( A-U ) “ ( 1 ••U) ) A 1640 

Z2s • 5* < A-U ) *AMXL2-AMXL I - • 5* ( I • -U ) *0MXL2+0MXL 1 A I 650 

130 CONTINUE A 1 660 

IF (A*LT#E) GO TO 140 A1670 

IF (ABS(A-l • I.LT.E) go TO 150 A1680 

140 H=-*5 A1690 

Qal,0 A1700 

Zl sU*AL 0G(U)-»5*U-*5* ( I •-Ul^OMXLl -♦••5»0MXL-.5 Al7l 0 

GO TO 160 A 1720 

150 HaO.O A1730 

0*H 

Zl s-OMXLl A 1750 

GO TO 160 A I 760 

160 TANTHO! U)«CLl*(Zl/( 1 • -Q* A ) - 1 • -ALOG ( U » -H ) /P 1 / ( A + 1 • )/2«0 A I 770 

IF (IC<Y.GT«1#AN0«L.LT»ICKY) GO TO 120 A 1 780 

IF CICKY«EQ*1) GO TO 100 A1790 

DO 170 J« 2 * ICKY AieoO 

170 TANTHOI 1 )sTANTHO( n-fTANTHO( J) A1810 

180 CONTINUE A I 820 

IF (ABS<CMB).LT#E 1 TANTHO»E A103O 

YP=10.»»10 A1840 

YPPalO.##lO A 1850 

YUPa- 1 /TANTHO A 1 86C 

YLPa- 1 /TANTHO A 1 670 
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790 • X f XU ( I > # YU ( n * XUC f YUC ♦ YUP • XL < I) 


770* XtYf YP*YPP.XC# YC 


GO 

TO 

240 

A2130 

GO 

TO 

250 

A2140 

GO 

TO 

260 

A2I50 

GO 

TO 

270 

A2160 

GO 

TO 

270 

A2170 


190 CONTINUE 
I ■ I 

IF ( ABS(CMB)«GT.0«l#*6> PRINT 
1 «YU( I) *XLC*YLCfYLP 
IF (A0S<CMB).LE*Otl*#6> PRINT 
XS.00025 
200 CONTINUE 

IF (PROFlLE»EQ»lOH4-DiGlT ) GO TO 210 
IF {PROFILE*EQ*10H4-D!GITMOD) GO TO 220 
PRINT 800 
GO TO 230 

210 Y»0#29690»SQRT (X )-0# 1 2600»X-0 • 35 1 60»X**2+0 • 284 30»X#<^3-0 • 10 15*X<H*4 
YP=•5^^•2969/SQRT(X)••l26-^•♦•35l6^^X+3•*•2e43*X^^X-4t I 015^^X^H^3 
YPPS-. 5#*5*. 2969/SORT (X»*3> -2.* #351 6+2# ♦a. *• 2043»X-3t »4*t • 1015»X»X 

GO TO 230 

220 YsA0*X»*«5+A I *X+A2»X**2+ A3»X«*3 

YP=«5»A0/X**#5+A 1 •f2.*A2*X+3*»A3#X##2 
YPP=-*25*A0/X*#1 •5+2«*A2+6«*A3*X 
230 CONTINUE 

Y»Y*T0C/«2 

VP=yP*T0C/#2 

YPP=YPP*T0C/*2 

IF (ABS(CMa>«LT«E) CM»0«5 

XC=X*CHD 

YC=Y»CHO 

IF (CAMBER«F0*10h2-0!GIT ) 

IF (CAM8ER.EQ«l0H3-01GIT ) 

IF (CAMBER.EQ. 10H3-DIGI TREF ) 

IF (CAMBER»EQ»lOH6-SERIES ) 

IF (CAMBER#E0*10H6A~SERIES ) 

PRINT 760 
GO TO 440 

240 YCMB=CM0* (2.0*CM*X-X*X)/CM»*2 
TANTH=2#*CMB*< I •-X/CM)/CM 

IF (X»GT»CM) YCMB=CM0*< I *-2« #CM+2 •»CM*X-X*X )/( I 
IF (X»GT#CM) TANTHs (2#*CM-2« *X )*CMB/ < 1 •-CM)*»2 
F=SQRT( 1 •+TANTH*»2) 

THP=-2#*CMB/CM#*2/F**2 

IF (X.GT.CM) THPa"2**CMB/( 1 .-CM)##2/F**2 
GO TO 440 

250 YCMB = CMB* < X**3-3* *CM*X»*2+CM»»2* ( 3« -CM ) #X ) /6 • 

TANTH = CMB* ( 3»*X^**2-6«‘»CM#X+CM»*2«’ ( 3* -CM ) )/6# 

IF (X«GT*CM) YCMB=CM0#CM**3» < 1 •-X)/6# 

IF (X*GT*CM) TANTHs-CMB*CM»t3/6« 

FaSQRT ( I •+TANTH*#2 ) 

THPsCMB* (X-CM)/F#»2 
IF <X#GT*CM) THPsO.O 
GO TO 440 

260 YCMBaCMB* ( ( X-CM > ♦♦3-K20K I ♦ ( I -CM ) ##3^^X-CM*#3*X-fCM*#3 >/6 
TANTHsCMB* ( 3»» (X-CM ) #»2-K20K 1 ♦ ( I -CM ) »»3-CM»»3 ) /6« 

IF <X*GT «CM ) YCMD = CMB* <K20K 1 » (X-CM ) »»3-K20Kl ♦ ( 1 -CM ) ♦*3*X-CM*«3*X+C 
lM*#3)/6 

IF (X.GT.CM ) TANTHsCMB* ( 3*K20Kl ♦ ( X-CM > *#2-K20K I # ( I -CM ) #*3-CM*#3 ) /6 
FsSQRT( 1 . + TANTH<^»2) 

THP=CMB# (X-CM)/F*#2 

IF (X.GT.CM) THPsK20Kl#CMB* (X-CM)/F*»2 
GO TO 440 
270 L=0 

A(l )*A3 
CLI'( I jaCLlS 
280 LsL+l 

Ad )sA(L) 

CLl < I )sCLI (L ) 

XCaX»CHD 
YC»Y»CHD 


Aiaao 

A I 890 
A I 900 
A19I0 
A1920 
A 1930 
A 1940 
A 1950 
A I960 
A 1 970 
A I 980 
A 1990 
A2000 
A2010 
A2020 
A2030 
A2040 
A2050 
A2060 
A2070 
A2060 
A 2 090 
A2100 
A21 10 
A2I20 
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A2160 

A2190 

A2200 

A2210 

A2220 

A2230 

A2240 

A2250 

A2260 

A2270 

A2260 

A2290 

A2300 

A23I0 

A2320 

A2330 

A2340 

A2350 

A2360 

A2370 

A 2 360 

A 2 390 

A2400 

A2410 

A2420 

A2430 

A2440 

A2450 

A2460 

A2470 

A2480 

A2490 

A2500 

A2510 

A2520 
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IF 

IF 

IF 

IF 


290 

300 

310 


XLLsK*AL0G(X ) 

Qs t «0 

IF (ABS(l.-A».UT.E.ANO»ABS<l.-X).Lr.E) 60 TO 330 
( A«LT.E»AND. ( 1 »-X »*LT»E ) GO TO 340 
(ABS(A-X)»i_T.E» GO TO 290 
( ABS( 1 .-Xj^LT.E ) GO TO 310 
(ABS(A-l. l.L-T.E) GO TO 320 
Va- ( A-X 1 /ABS ( A-X ) 

OMXL« ( I ••X )*AL0G( I *~X ) 

AMXU= ( A-X ) tALOG ( ABS ( A-X ) ) 

OMxLls-ALOG(U-X)-l • 

AMXLI=-AL0G(ABS<A-X) )+V 
0MxL2*l i/da-X) 

Amxl 2*1 t/CA-X) ^ 

Zs«5*( A-X)*AMXL-#5* < l • -X ) ♦0MXL-*25# ( A-X ) ♦♦2-f# 25* < 1 *-X)*#2 
Zla*5*( (A-X)*AMXLI-AMXL- ( I *-X )*OMXLl+OMXL+ ( A-X )- ( 1 •-X) ) 

Z2**5* < A-X )*AMXL2-AMXL 1 -#5* ( 1 *-X > *0MXL2+0MxL I 
IF (A.LE.E) GO TO 300 

Ga-( A*A* ( •5*AL0G( A )-0#25 )+0«25>/ < I .-A ) 

Ha (0*B» n #-A )**2*AL0G< 1 »-A )-0«25»( I t-A )*#2 )/ ( I a-A ) +G 
GO TO 350 

2s-.5»< I •-X)<H^2*ALOG( I •-X)‘»-0*25*< 1 *-X)**2 

21 a-*5<M I «-X 1* ( -ALOG ( 1 • -X > - 1 • ) + • 5* ( I t-X )*AlOG( I #-X )-r*5* ( I *-X ) 
Z2=-ALOGU •-X)-0#5 

Ga-( A»»2 * ( .B^AlOG (A )-0#25 )+0»251/ < I • -A ) 

Ha (0i5<M 1 .-A )**2*AL0G( I *-A ) -0.25* U • -A ) **2 ) / ( 1 t-Al+G 
GO TO 350 
Ga-*25 
Ha-*B 
GO TO 350 
CONTINUE 

Ga-(A*»2*(#5*ALOG(A)-0t25)+0«25)/< I »-A) 

Ha (0«5» < 1 *-A )*»2*AL0G( 1 #-A )-0*25» ( 1 #-A )»*2 )/( I •-A)+C 
Za#5*(A-l • )##2»ALOG( ABS(A-l • ) >-0.25*(A-l • )**2 
21 a- (A- I • )#ALOG ( ABS ( A-l • ) ) 

22«-l0f <H^I0 

GO TO 350 
320 G*0#0 
H*G 
QaG 

2a-(l.-X)*ALOG<U-X) 

21 aALOG( 1 *-X ) + l • 

Z2*-l •/< I .-X ) 

GO TO 350 
330 2*0.0 
G*Z 
H*Z 
Q*Z 

21a-lO.«#lO 
22*-l0.#*10 
GO TO 350 
340 Ga-.25 
Ha-,5 
0*1.0 
Z*-.25 
21*0.0 
22*-l0.#*t0 
GO TO 350 

350 YCMB(L)«CLl*<Z/< I •-Q^^A)- 
XSV»X 

IF (X.LT.0*005) X*0.005 

TANTH(L)*CL|»(Z1/(1#-‘3*A)-1 .-AL0G(X 1-H)/PI/( A+l . )/2.0 
XaXSV 

IF ( IF6XA(L).EQ*l ) TANTH(L)«-5. 

IF (X«GT.O«005) GO TO 360 
YCPr (L 1 *0.0 
GO TO 380 


-XLL4*G-H#X)/PI/(A+1 • )/2. 


A2530 

A2540 

A2550 

A2560 

A2570 

A2580 

A2590 

A2600 

A2610 

A2620 

A2630 

A2640 

A2650 

A2660 

A2670 

A2680 

A 2690 

A2700 

A27I0 

A2720 

A2730 

A2740 

A2750 

A 2760 

A2770 

A2780 

A2790 

A2600 

A2810 

A2B20 

A2830 

A 2840 

A2850 

A2860 

A2870 

A2860 

A2890 

A2900 

A29I0 

A2920 

A2930 

A2940 

A2950 

A2960 

A2970 

A2980 

A 2990 

A3000 

A3010 

A3020 

A3030 

A3040 

A3050 

A3060 

A3070 

A3060 

A3090 

A3I00 

A31 10 

A3120 

A3130 

A3140 

A3150 

A3I60 

A3170 

A3IB0 

A3190 
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360 CONTINUE A3200 

IF ( ABSU .-X ).GT#E) GO TO 370 A32l? 

YCP2(L )«l »/E A3220 

GO to 380 A3230 

370 P|AaPl#(A+l * )*2« A3240 

YCP2 (L)aCLl* <Z2/( I t-Q^^A )-l »/X )/P I A A3250 

380 CONTINUE A3260 

MODIFIED camberline OPTION A3270 

IF <CAMBER.EQ«K6AS) GO TO 390 A3280 

GO TO 410 A3290 

390 YCMBCL)=YCMB(L)*0.97948 A3300 

TANTH CL )sTANTH(L > #0.97948 A3310 

YCP2CL >aYCP2<L)#0.97948 A3320 

IF ( ABS(A-.8).LT.E) GO TO 400 A3330 

PRINT 780 A3340 

read 600* NPWIPE A3350 

IF (K0N.EQ#3) KON=0 A3360 

GO TO 20 A3370 

400 CONTINUE A3380 

IF (TANTH<L>.LE.“#24521#CLI ) YCMBIL > *0 • 2452 1 *CL I # ( I .-X ) A3390 

IF (TANTH(L)«LEo-#2452l#CLl ) YCP2(L)*0.0 A34 q0 

IF (TANTH(L).LE«-*2452l#CLl ) T ANTH ( L) =-0 • 2452 I #CL I A3410 

IF (TANTH(L).LE«-#24521#CLI ) IFeXACL)*! A3420 

410 CONTINUE A3430 

IF ( ICKY.GT* I .AND.L.LT# ICKY ) GO TO 280 A3440 

IF (ICKY.EQ.l) GO TO 430 A3450 

DO 420 J=2*ICKY A3460 

YCMBC 1 )=YCMB<1 )+YCM8(J) A3470 

TANTHCl )=TANTH(1 )+TANTH(J) A3480 

,YCP2 ( I ) =YCP2 ( I >+YCP2 ( J ) A3490 

420 CONTINUE A3500 

430 CONTINUE A3510 

F=SQRT( 1 •+TANTH##2) A3520 

THPsYCP2/F**2 A3530 

440 CONTINUE A3540 

IF (X.GT.XM) GO TO 550 A3550 

IF (ABS(X-XM)#LT*E) GO TO 550 A3560 

SINTHsTANTH/F A3570 

COSTH^l */F A3580 

I a 14-1 A3590 

XU( 1 )sX“Y#SI NTH A3600 

YUC I )sYCMB+Y#COSTH A361 0 

XLC I )«X+Y#SINTH A3620 

YL< I )=YCMB-Y#C0STH A3630 

IF (ABS(X-XA(K) ).GT.#l##6) GO tO 450 A3640 

XAUCK )sXU( I ) A3650 

YAUCK )«YU( I ) A3660 

XAL(K)=XL( I ) A3670 

YAL(K)=YL( I ) A3680 

KsK+1 A3690 

450 CONTINUE A3700 

XUCaXUC I )#CHD A3710 

YUCsYUC I )#CHD A3720 

XLC*XL( 1 )»CHD A3730 

YLC»YL( I )#CHD A3740 

IF (ABS(CMB).LE*0.1##6) GO TO 460 A3750 

YUPeO.O A3760 

VLPsYUP A3770 

IF (ABS(TANTM).LT*0.1##10) GO TO 460 A3780 

YUP= < TANTM#F + YP-TANTH*Y#THP )/<F-YP#TANTH-Y#THP) A3790 

YLPs (TANTH#F-YP+TANTH»Y#THP )/ (F+YP#TANTH4^Y#THP) A3800 

460 CONTINUE A3810 

IF <X*LE»0.0975 ) FRAC=0.25 A302O 

IF (X.LE.0.00225) FRAC*0.025 A3830 

IF (AB5<CMB)*GT.0#l##6) PRINT 790* X • XU ( I ) • YU < I ) • XUC * YUC • YUP • XL ( I ) A3840 


o o o 
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I f YL( n *XLCtYLCfYUP ^ ^ 

IF UHS(C,Mfci)tLE«0#l*#6) PRINT 770. X . Y . YP . YPP . XC • YC 
Xax->F RftC^DX 
FRACal #0 

IF (AB5(X-xM)«UT.E) go to 470 
IF (X*LT»XM) GO TO 200 

PROFILE - X GE XM 

X»XM 

470 CONTINUE 

IF (PROFiLEtEQi 10H4-OIGI T ) GO TO 4B0 
IF (PROFlLE*.iQ*10H4-*DlGITWOO) GO TO 490 
PRINT 800 
GO TO 500 

480 Y»0«29690»SQRT f X >-0* 1 2600»X-0 .35 1 60*X#*2 + 0 •28430#X**3-0 • 1 0l5<fX##4 

YPa.3#«2969/SQRT (X)-* l26-2#*«35l6#X+3.*.2a43»X*X-4«». 1015»X**3 

YPP = ~,5#*5#*2969/SQRT ( X^**3 > -2« * • 351 6 + 2 • * •2B43*X-3#»4 • ♦ ♦ 1 01 5*X*X 

GO TO 500 

490 Y = D0+DI * < 1 #-X )+D2* ( I #-X ) **Z + 0:^* < I «-X >**3 

YPa-Dl -2**D2* (1 .-X ) -3**D3» < 1 *-X ) 

YPP=2**^D2+6«*D3* ( 1 *-X ) 


A3B50 
A 3860 
A3870 
A38B0 
A3B90 
A39O0 
A391 0 
A3920 
A3930 
A3940 
A3950 
A3960 
A3970 
A3980 
A3990 
A4000 
A4010 
A4020 
A4030 
A4040 
A4050 
A4060 


IF 

IF 

IF 

IF 


GO 

TO 

510 

A4130 

GO 

TO 

520 

A4140 

GO 

TO 

530 

A4150 

GO 

TO 

540 

A4160 

GO 

TO 

540 

A4170 


500 CONTINUE 

Y=Y*T0C/*2 
YPrYP*T0C/*2 
YPP=yPP*T0C/.2 
XC=X*CHD 
YC=Y*CHD 

IF (CAM0ER»EQ*IOH2-DIGIT > 

(CAMBER.EQ* IOH3-D IGIT ) 

( C AMBER •EQ# I 0H3-DIGI TREF > 

(CAMBER.EQ* 10H6-SERIES ) 

( CAMBER .EQ# I 0H6A-SERIES > 

PRINT 760 
GO TO 560 

510 YCMB=CMB* (2«0*CM#x-X*X )/CM**2 
TANTH=2*»CMB» (1 .-X/CMl/CM 

IF (X.GT.CM) YCMB=CMB*( I • -2 • #CM+2* *CM^^X-X#X )/( I •-CM>##2 
IF “(X*GT#CM) TANTH=(2**CM-2.*X)*CMB/n 
FsSQRT < I • + TANTH*#2 ) 

THP=-2**CMB/CM#*2/F**2 

IF (X.GT.CM) THP*-2#^fCMB/( I • -CM ) #»2/F**2 
GO TO 560 

520 YCMBaCMB* (X<f»3-3##CM*X**2+CM#*2# (3*-CM)#X )/6# 

TANTH = CMb*< 3#*X**2-6*»CM#X + CM*»2*(3*-CM > )/6# 

IF <X#GT*CM) YCMBsCMB#CM»#3* ( 1 *-X)/6* 

IF (X.GT.CM) TANTH=-CMB*CM»»3/6. 

FsSQRT < I • + TANTH**2 ) 

THPsCMB^t (X-CM)/F*»2 
IF (X*GT#CM) THPsO.O 
GO TO 560 

530 YCMB«CMB*( (X-CM)#*3-K20Kl*( 1 -CM ) #*3*X-CM*#3* X+CM## 3 ) /6 
TANTH = CMB» ^ 3 ^* <X-CM) #*2-K20< 1 ♦ U -CM ) #*3-CM**3 ) /6«« 

IF <X*GT.CM ) YCMa = CMB*(K20Kl » ( X-CM ) » »3-K20< 1 *( I -CM ) 4»3*X-CM**3*X + C 
lM#«3)/6 

IF ( X»GT*CM ) TANTHsCMB* ( 3»K20K1 ♦ ( X-CM ) #«-2-K20K I * ( 1 -CM )*#3-CM#*3 )/6 
FsSQRTI l. + TANTH♦^^2) 

THP = CM0^MX-CM)/F*»2 

IF (X*GT#CM ) THPSK20K1 <^CMB» (X-CM )/F**2 
GO TO 560 
540 GO TO 270 
550 CONTINUE 
560 CONTINUE 

S|NTH«TANTH/F 
COSTH*l ./F 
IM + l 


A4070 
A4080 
A4090 
A4100 
A41 10 
A4120 


A4180 

A4190 

A4200 

A42IO 

A4220 

A4230 

A4240 

A4250 

A4260 

A4270 

A4280 

A4290 

A4300 

A4310 

A4320 

A4330 

A4340 

A4350 

A4360 

A 4 370 

A43B0 

A4390 

A4400 

A4410 

A4420 

A4430 

A4440 

A4450 

A4460 

A4470 

A4480 

A4490 

A4500 
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XU( I )»X-Y*SINTH A4510 

YU( I »=YCMB+Y*COSTH A4520 

XL ( t >*X + Y*S1NTH A4630 

YL ( I >=YCMB-Y*COSTH A4540 

IF (ABSIX-XA (K ) )tGT.. 1**6) GO TO 570 A4550 

XAIXK ) aXU( I ) A4560 

YAU(K)=YU( I ) A4570 

XAL(K)=XL( I > A4580 

YAL(K)=YL(I) A4590 

K»K+1 A4600 

570 CONTINUE A4610 

XUCsXUC I )*CHD A4620 

YUC»YU< I )*CHD A4630 

XLC»XL( I )*CHD A4640 

YlC«YL ( I )*CHD A4650 

IF (ABS(tMB).LE*0.I**6) GO TO 580 A4660 

YUPsO.O A4670 

YlPsVUP A460O 

IF (ABS(TANTH).UT.0«l**10> GO TO 580 A4690 

YUPsTANTH* (F+YP/TANTH-Y*THP )/(F-YP#TANTH-Y#THP ) A4700 

YLP=TANTH» (F-YP/TANTH + Y#THP )/ (F+YP*TANTH + Y»THP ) A47l 0 

580 CONTINUE A4720 

IF (ABS (CMB ) •GT*C#l **6) PRINT 790« X • YU < I ) * YU ( 1 ) • XUC % YUC • VUP • XL < I ) A4730 

1 #YL( 1 ) *XLC%YLC%YLP A4740 

IF (ABS(CMLi).LE*0*l^*6) PRINT 770% X % Y % YP % YPP % XC % YC A4750 

X=x+DX A4760 

IF (X«LE.l#0) GO TO 470 A4770 

PUNCH 600% (TITLEU ) %1 =1 %3) A470O 

PUNCH 810% ( XAu ( I ) % I =1 % 32 ) A4790 

PUNCH 810% (YAU( I )% I -I %32) A4791 

PUNCH 810% (XAL( I )% 1 = 1 %32) A4792 

PUNCH 810% (YAL< I )% 1*1 %32) A4793 

CALL PLOT (XU%XL% YU%YL% n A40OO 

GO TO 20 A4810 

A4820 

590 FORMAT ( I 3/ ( 8F 1 0 • O ) > A4630 

600 FORMAT ( 8A 1 0 ) A4B40 

610 FORMAT (IH %8AlO> A4850 

620 FORMAT (lOHlPROFILE %A10%10H CAMBER %A10> A4S60 

630 FORMAT (7FlO%0%AlO) A4870 

640 FORMAT (8FlO%0) A4880 

650 FORMAT ( / I OX % A I 0 % 1 OX % A I 0/ ) A4890 

660 FORMAT (19H PROFILE PARAMETERS/5H T/C* %F I 0»5/ 1 2H L %E %RAD I US= % F 1 0 • 5 A4900 

1/18H BASIC X INTERVAL*%F10%5/7H CHORD* %F I 0% 5/ ) A49lO 

670 FORMAT <19H PROFILE PARAMETERS/5H* T/C* % F I 0% 5/ I 2H L %E %RAD I US= %F 1 0 *5 A4920 

1/18H BASIC X 1NTERVAL*%F10%5/7H CH0RD= % F 1 0% 5/35H POSITION OF MAXIM A4930 
2UM THICKNESS% XMs %F I 0# 5/ 1 3H CONSTANT D1=%F10%5/) A4940 

680 format (23H CAMBER LINE PARAMETERS/ 1 6H CAMBER ( YCMAX ) s%riO*5/28H P A4950 
lOSlTlON OF MAXIMUM C AMBER* % F 1 0 • 5/ ) A4960 

690 FORMAT (23H CAMBER LINE PARAMETERS/2 1 H CAMBER PARAMETER Kl = %F10%5/ A4970 

140H POSITION OF ZERO CAMBER LINE CURVATURE =% F 1 0 • 5/ ) A4980 

700 FORMAT (23H CAMBER LINE PAR AMETERS/2 1 H CAMBER PARAMETER Kl*%riO%5/ A4990 
140H POSITION OF ZERO CAMBER LINE CURVATURE* ♦ F I 0 • 5/6 1 H RATIO OF AFT A5000 
2 TO FORWARD CAMBER LINE CURVATURE FACTQR% K20X 1 * % F I 0 %5/ ) A5010 

710 FORMAT (23H CAMBER LINE PARAMETE RS/ 7X % 3HCL I % 9X % 1 H A ) A5020 

720 FORMAT (2FI0.3) A5030 

730 FORMAT ( 1 OH AO % 1 % 2 % 3 * % 4F 1 3% 6/ 1 OH DO % I % 2 % 3a % 4F I 3* 6/4H RC*%F13%3//) A5040 

740 FORMAT (9X%3HX/C% 10X%3HY/C%3X%5HDY/DX%6X%7HD2Y/DX2%22X% 1HX% 12X% IHY A5050 
l/j A5060 

750 FORMAT ( I 1 6H0UNC AMBEREO UPPER SURFACE VALU A5070 

LOWER SURFACE VALUES /5X% A5080 

23HX/C1 7X%4HXU/C6X%4HYU/C5X% 7H XU 3X % 7h YU 3X % 7HDYU/DXU I 3X % ^H A5090 
3XL/C4X%4HYL/C5X%7H XL 3X % 7H YL 3X % 7H0VL/ DXL > A5I00 

760 FORMAT (35H BAD HOLLERITH CAMBER SPECIFICATION) A51I0 

770 FORMAT ( 4F 1 3 • 6 % I OX % 2F 1 3 »6 ) A5I20 

700 format (52H MODIFIED CAM3ER LINE OPTION ONLY ALLOWED IF A=0%8 > A5130 

PAp.13 13 

OF 


1 


M' 




APPENDIX 

790 FORMAT < F I 0 . ‘ij • I OX • ^4P I 0 • 5 • E I I «2 • 6X t4F I 0 . 5 • E I I • 2 ) 
aOO FORMAT ( 36H UAO HOLLERITH PROFILE SPECIFICATION) 

BIO FORMAT (aFlOi^i) 

END 

SUBROUTINE PLOT ( XU ^ XL t YU • VL • I ) 

COMMON /MAIN/ YSTART ( 3 ) f CHDtK 

DIMENSION XU(l). XL ( 1) . YU(l)t YL ( 1 ) i X(460)* Y(4b0) 

DIMENSION TITLEMB)* TITLE2(B) 

READ 30f (TITLEl (N) •N-l#a) 

IF (M0D(K»3).EQ*1 ) CALL CALPLT (l#0%0*0»-3) 

HGT = Ot 14 
L* I 

DO 10 Nal • I 
X(N)«XU(N) 

Y(N)aYU(N) 

X(J+N)aXL(L) 

Y ( I+N) »YL (L ) 

10 L*L-I 
Ms2*l 
XPG*lO#0 

XX = XPG/2«0 - 1 #5* (6 */7 •■*HGT ) 

XOVaO.O 
XT1C=1 #0 
YPGs2*0 
YDV=0*0 
YTIC«l «0 
X<M+1 )=0#0 
YIM^I )=-0.1 
X(M4-2 ) = 1 *0/XPG 

CALL^ AXES^^ 0 ! , YST ART ( K ) 1 90 . • YPG» Y ( M+ 1 ) • Y ( M+2 ) • YT I C • YDV . I H • HGT . I ) 
CALL AXES (0 t » YSTART (K > «0. *XPG«X (M+l ) iX CM+2 ) *XT ICfXDVf I H tHGT*-!) 
YlABEL=YSTART (K)- 2#b*HGT 

call NOTATE ( XX ♦ YLABEL • HGT • 3HX/C • 0 • « 3 ) 
ylabel=ylabel-i *5»HGT 

CALL NOTATE ( 0 • 0 « YL ABEL • HGT ♦ T I TLE I < • 0 • 90 ) 

YS=YSTART (K )+l .0 

NOTATE ( — * 92 • Y 5 • HGT ♦ 3HZ/L • 0 • 0 • 3 ) 

CALL CALPLT ( 0 • 0 , YS T ART ( K ) • - 3 ) 


LAP = 0 

call linplt 

CALL CALPLT 
IF (K*LT*3) 
KaO 

CALL NFWAME 
20 CONTINUE 

return 


<X f YfM« 1 «LAP«0» I «0 ) 
(O.Ot-YSTAHT <K ) f-3) 
GO TO 20 


A5I40 
A9160 
AB160 
A5I70- 
0 10 


B 

a 

B 

B 

B 

B 

B 

B 

B 


30 FORMAT 
END 


(8A10) 


20 
30 
40 
50 
60 
70 
BO 
90 
100 
B no 

B 120 
B 130 
B 140 
B 150 
B 160 
B ItO 
B Ido 
B 190 
B 200 
B 210 
B 220 
B 230 
B 240 
8 250 
B 260 
B 270 
B 260 
B 290 
B 300 
B 310 
B 320 
B 330 
B 340 
B 350 
B 360 
B 370 
B 380 
B 390 
B 400 
B 410 
B 420 
B 430 
B 440 
B 450 
B 460- 
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TABLE n.- SAMPLE COMPUTER PRINTOUT OF ORDINATES 
FOR SYMMETRIC AIRFOIL 

PRnriM' 2-nroiT 


ppnr II *- PARAMFTPP'’. 

T/r= . lanoo 
L.r •pa'U!js= -o.oojoo 
BASIC X INTF»VAL^ .OlOOO 
CHORHs A.000J3 

I- INC PAPAMF Tri>s 
CAMBFRIYCMAXI - 0.00000 


TION nr maximum. CAMH^Rs 

0.00000 




x/c 

Y/C 

OY/OX 

n2Y/nx2 

x 

Y 

o.ooooon 

0.000090*00000.000000*00000 .000000 

0.000000 

0.000300 

. J00?50 

. 002 700 

4 . 557476 - 

11266.954512 

.OJ15JO 

.015786 

.000500 

.003045 

3.007521 

-3953.752904 

.0030UU 

.023673 

,000750 

.004R2? 

3.176460 

-2168.672154 

.004500 

.023930 

. 001000 

.005557 

2.740619 

-1408.741253 

• jJ6 j jO 

.033345 

.001250 

.005203 

2.443153 

-1008.132658 

.007500 

.037220 

.001500 

.006755 

2.223546 

-767.01259) 

.009000 

.0407 1 3 

.001750 

.007319 

2.052843 

-605.757344 

.010500 

.043915 

.002000 

.007514 

1.914224 

-495.336313 

.012000 

.045883 

.002250 

.005279 

1.501214 

-41 7.699726 

.013500 

.049673 

.002500 

.00571*^ 

1 . 70474 v5 

-356 .699366 

.015000 

.052300 

. j 05 oOO 

.012213 

1. 181943 

-126.380823 

.030000 

.073279 

,007500 

.014549 

.949756 

-68.980403 

.045000 

.089391 

.OlOOUO 

.017037 

.510932 

-44.946745 

.060000 

. 102222 

.012500 

.01593^ 

.7158*»2 

-32.275893 

.075000 

.113634 

.015000 

.020407 

.^44439 

-24.6485U 

.090000 

.123320 

.01 7500 

.022179 

.4Q0478 

-19 .641545 

.105000 

.133J74 

. 020000 

.023590 

.545954 

-16,147243 

.12000J 

. 141587 

.022500 

.024915 

.505963 

-13.594517 

.135000 

.149491 

.0?5u00 

.026147 

.477496 

-ll .663353 

•150000 

. l*=6083 

.027500 

.027306 

.450291 

-10,160005 

.165000 

. 163837 

.050000 

. u254J 1 

.426442 

-8.962635 

.130000 

.170439 

.032500 

.029441 

.404291 

-7.990525 

.195000 

. 176644 

.035000 

.030430 

.35045 

-7.155410 

.210000 

.182578 

.037500 

.031374 

.369240 

-6.517307 

.225000 

.188243 

.340000 

.032277 

.353676 

-5.949025 

.240000 

.193663 

.0^2500 

.003143 

.339427 

-5,462718 

.255000 

.198360 

.045000 

.033974 

. 326 30 5 

-5.042677 

.270000 

.203852 

.047500 

.03477F 

. 314 165 

-4,6 76 8 61 

.235000 

.203654 

.0*^0000 

.035*^47 

.302556 

-4.3^5004 

.300J00 

.213281 

.052«;00 

.036291 

.292 350 

-4.072424 

.315000 

.217745 

.055000 

.037009 

.282495 

-3.820528 

•330000 

.222055 

.057500 

.037704 

.273233 

-3.595465 

.345000 

.226222 

.040000 

.035-'?6 

.264501 

-3.393365 

.360000 

.230255 

.0F«2500 

.039027 

.256250 

-3.211047 

.375000 

.234160 

.045000 

.039447 

.248432 

-3.045876 

.390000 

.237945 

.067^00 

.040269 

.241008 

-2.895651 

.405000 

.241615 

.070000 

.040543 

.233943 

-2.758521 

.420000 

.245177 

.072500 

.041439 

.227205 

-2.632922 

.435000 

.243635 

.375000 

.'j4199<' 

.220769 

-2.517521 

.450000 

.251994 

.077500 

.042543 

.214610 

-2.411176 

.455000 

.255259 

.050000 

.043072 

.205707 

-2.312906 

.483000 

.2*^8434 

.0fl?‘'00 

.043557 

.203040 

-2.221863 

.495000 

.261522 

.055000 

.044055 

. 197592 

-2.137308 

.510000 

.264526 

.057500 

.044575 

. 192 345, 

-2.058598 

.525000 

.267* 50 

.ooonoo 

.045050 

.157295 

-1.985171 

•540000 

.270297 

.092500 

.045512 

.182419 

-1 .916530 

.555000 

.273070 

.095000 

.045962 

.177705 

-1 .552239 

.570000 

.275771 

.097^00 

.046400 

.173144 

-1.791910 

.58^000 

.278402 

.100000 

.046525 

.165746 

-1 .735200 

.600000 

.280966 


''.RIGINALFAG®® 
OF POOR QUALiriil 


TABLE n.- Continued 


x/c 

.noooo 

.i?oooo 
, I ^0000 
.UOOQO 

.isoooo 

. lf>0000 
.170000 
.IBOOOO 
.100000 
.?00000 
.?10000 
.220000 
.2*^0000 
. ?<»0000 
.7*^0000 
.760000 
.770000 
.7B0000 

. -^ooooo 

.300000 
.310000 
,320000 
,330000 
.3A0000 
.350000 
.360000 
.370000 
.330000 
. 390000 
.AOOOOO 
.AlOOOO 
.470000 
.430000 
. 440000 
, 4 3000 
,4^3000 
,470000 
,4.U)000 
,4^0000 
OOOOO 
.MOOOO 
,‘=70000 
. 5300 JO 
.*'43000 
, *^‘•,0000 
.*■•60003 
• 5 7'Ju UO 
.‘■BO 000 
.*'*■‘0000 
. A.OOOOO 
13303 
>0300 
.43 1003 
.640000 
.'•rooOO 
^ 3003 
.470000 
A0300 
.'03300 

. 000 on 

.710303 
, V0003 

, ^3 00 on 

, 740000 
. ■'**'0 300 
. 760000 
. '*70 0 00 
, '>33003 


Y/C 

.048437 
.04983? 
.3*=ll93 
.057380 
.053452 
. J5441P 
.055787 
.0'’6066 
.0*=6760 
.057375 
.057916 
,0‘=R386 
.058790 
.059131 
.059412 
.059637 
.059807 
.059926 

. OC999C 

.060017 

.059904 

,059978 

.0*'98?0 

.05967? 

.059486 

.059763 

. 0*^9005 

.0*=87l? 

,058387 

.0^8330 

.057643 

.05777c 

. 0*^6780 

, 0*=6307 

. 80 7 

.:)55787 

.'^54737 

.053560 
. 3 *= 7 94 3 
. >‘-?>9R 
.351635 
.050951 
,350’>48 
.049574 
.34878? 
.348321 
.047743 
.046447 

.045634 
,0 44 834 
.u43 r,8 
,343096 
.04771 8 

.041375 

.34041 7 
.1)3949 4 
.030557 
.33 7*0'^ 
.336539 
.33«=4C9 

. 3344' (. 

. 3^3658 

.932637 

.331603 

.979494 

.378473 


DY/DX 

. 15241 3 
. 137841 
. 124699 
. 1 17742 
.101782 
.091671 
.087293 
.O735C4 
.065379 
.057703 
.050475 
.043650 
.037189 
,031061 
.025238 
.019694 
.014408 
,009352 
,004538 
-.000078 
-.004500 
-.008740 
-.012808 
-.016716 
-.020472 

-.024085 

-.027562 

-.03091 I 

-.034137 

-.037748 

-.040249 

-.043145 

-.045942 

-.348645 

-.051257 

-.053734 

-.056278 

-.058595 

-.0608*’8 

-.06311 I 

-.06526 6 

-.067358 

-.069389 

-.07137>3 

-.073782 

-.075149 

-.076967 

-.078739 

-.080467 

-.0821*^4 

-.083803 

-.085415 

-.086994 

-.088541 

-.0900*^ 9 

-.091550 

-.093016 

-.394459 

-,09»*Rn2 

-.097787 

- ,398676 

-. 10JJ50 

-.10141 1 

-.107763 

-.104106 

-.105442 

-.106775 

- . 108104 


D2Y/DX2 

-1,538890 

-1 .30097? 

-1,251356 

-1.143134 

-1.051434 
-.972731 
-.904421 
-.844538 
-.791579 
-.744372 
-.701995 
-.663711 
-.628926 
-.597157 
-.568005 
-.541142 
-.516291 
-.493223 
-.471741 
-.451679 
-.432894 
-.415264 
-.398681 
-.383055 
-.368305 
-.354359 
-.341158 
-.328644 
-.3 16772 
-.305496 
-.294780 
-.204588 
-.274891 
-.265660 
-.256877 
-.248503 
-.240533 
-.2 32944 
-.225720 
-.2 18844 
-.2 12 304 
-.206086 
-.200179 
-.194573 
-.189257 
184272 
-.1 79 4 61 
-.174966 
-.1 70729 
-.166744 
-.163005 
-.15950^ 
-.156244 
-.153211 
-.150404 
-.147818 
-.145451 
-.143297 
-.1413*^4 
-.139618 
-.138087 
-.1 36757 
-.135626 
-.1 346 9? 
-.133951 
-.133403 
-.1 33044 

< -.137873 


X 

.660)00 
,720000 
.780000 
.840000 
.900000 
,960Ci0iJ 
1.020000 
1.080000 
1,140000 
1.200000 
I .260000 
1.320000 
1 ,380000 
1 ,440000 
1.500000 
I ,560000 
1.620000 
1.630000 
1.740000 
1,800000 
1 .860000 
1.920000 
1.980000 
2.040000 
2.100000 
2.160000 
2.220000 
2.280000 
2,340000 
2.400000 
2.460000 
2.520000 
2.580000 
2.640000 
2.700000 
2.760000 
2.820000 
2.880000 
2.940000 
3 . 0)0000 
3.060000 
3.120000 
3.180000 
3.240000 
3.300000 
3.360000 
3.420000 
3.480000 
3,540000 
3.600000 
3.660000 
3. 720000 
3.780000 
3.840000 
3.900000 
3.960000 
4.023000 
4.080000 
4.140000 
4.200000 
4.260000 
4.320000 
4,300000 
4.440000 
4.500000 
4.*'60000 
4.620000 
4.680000 


Y 

.290591 

.299291 

. 307161 

.314278 

.320710 

.326509 

.331725 

.336397 

.343562 

.344253 

.347496 

,350318 

.3*27 41 

,354787 

,356475 

.357821 

,3''0a43 

.359555 

.359971 

.360104 

. 359965 

.359567 

,358920 

. 358333 

.356917 

.355580 

.354J30 

.352275 

.3*^3323 

. 343181 

,34‘=85c; 

.343353 

,340680 

.337842 

.334844 

.331693 

.328392 

.324947 

.321362 

.317642 

.313790 

.309811 

.305708 

.301485 

.297146 

.292692 

.2G8129 

.283457 

,273601 

.273802 

.263323 

.263747 

.258574 

.2*’3308 

.247953 

.242501 

.236964 

.231340 

.275630 

.219334 

.213955 

.207994 

.201950 

.195824 

, 1896 18 

. 183337 

.176965 

.17)519 


TABLE li.- Concluded 


x/c 

Y/C 

DY/DX 

D2Y/DX2 

X 

Y 

7‘iOOOO 

.02733;> 

109433 

-.132000 

4.760000 

.163993 

nOOOOQ 

• 0?6’31 

-.110762 

-.133008 

4.800060 

.157387 

ftnooo 

.0?*>117 

-.112095 

-.133470 

4.860000 

.150701 

n ? 0 C )00 

.073980 

-.1 13432 

-.1 34033 

4.92J000 

.143936 

<^^0000 


-.114776 

-.134776 

4.900000 

. 137J89 

flAOOOO 

.021694 

-.116128 

-.135697 

5.040000 

.130162 

R50000 

.020826 

-.117491 

-.136794 

5.100000 

.123154 


.0193AA 

-.110068 

-.138068 

8.160000 

.115063 

870000 

.018K8 

-.120253 

-.139516 

5.220000 

. 10389U 

880000 

.016939 

-.121656 

-.141137 

5.280000 

. 101633 

800000 

.01 671*= 

-• 123076 

-.142931 

8.340000 

.094291 

oooooo 

.014477 

-.124515 

-.144896 

5.400000 

.086863 

910000 

.013228 

-.125974 

-.147031 

5.460000 

.079348 

*^?0000 

.011080 

-.127456 

-.149336 

5.520000 

.071746 

930000 

.010676 

-.120962 

-.151809 

5.580000 

.064053 

9<r0000 

.00937B 

-.130493 

-.154450 

5.640000 

.056270 

9*^0000 

.000066 

-.132081 

-.157258 

5.700000 

.048394 

960000 

.006737 

-.133639 

-.160232 

5.760000 

•040423 

O70000 

.005393 

-.135256 

-.163371 

5.820000 

.032356 

983000 

.004032 

-.13t90 7 

-.166675 

5.880000 

.024192 

990000 

.002684 

-.138591 

-.1 70143 

5.940000 

.015927 

000000 

.001260 

-.140310 

-.173775 

6.000000 

.007560 


TABLE in.- SAMPLE COMPUTER PRINTOUT OF ORDINATES FOR CAMBERED AIRFOILS 
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Ficure 2.- Variations in thickness-chord ratio for symmetric 4-digit airfoils 





NPCfl 2412 

- Variations in amount of camber for 12-percent-thick 4-digit airfoils with 2-digit camber line 













Variations in location of maximum camber for 12-percent-thick 4-digit modified airfoils 

with 3-digit camber line. 






(16-006) 

Variations in thickness -chord ratio for symmetric 16-series airfoils 




NfiCfi 16-212 

Variations in amount of camber for 12 -percent -thick 16-series airfoils with 6A-series camber line, 





Figure 11.- 2-digit and 3-digit camber lines 




